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Title:  ^ Multiping  Data  Processor  for  Target  Detection  * 

J * iT^  f . .. 

During  the  previous  quarter,  the  multiping  processor  design  was  finalized  nnf*  param- 
eters selected  for  the  detection  and  tracking  portions  of  the  multiping  processing  algorithm . 
Work  also  progressed  and  ip  continuing  on  obtaining  quantitative  performance  results  on 
t’.n  muhiping  processor  through  Monte-Carlo  techniques.-  Preliminary  performance  results 
will  be  obtained  by  January  30,  1970.  Figure  1 shows  the  l^lock  diagram  of  the  m-iCping 
receiver  configured  with  an  output  display  and  automatic  detection  device.*  An  analysis  of 
average  total  costs  for  multiping  sonar  decisions  was  made  and  will  be  useful  in  comparing 
com|Hiting  decision  algorithms.  This  analysis  is  presented, in  the  attached  Technical  Memo 
No.  1.  Technical  Memo  No.  2 ^evelops^he  Kalman  Alter  equations  and  doeertbes  the 
adaptation  oflhe  Kalman  filter  to  our  specific  track  detection  problem. 

During  the  next  quarter  quantitative  performance  results  of  the  multiping  processor 
will  be  obtained  via  Monte-Carle  methods  using  the  digital  tapes  generated  containing  noise 
and  signal  plus  noise  with  time  and  frequency  spreading  of  arriving  target  associated  paths . 
This  data  was  generated  using  the  ta pget/envi romnent  model  described  !a  the  Second  Pro- 
gress Report  dated  October  15,  1969. 
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Figure  1.  Multiping  Receiver  Bloch  Diagram 


Technical  Memo  No.  1 


January  2,  1970 
T.  G.  Kincaid 

Contract  No. : N00024-69-C-1268 


Title:  Average  Total  Costs  for  Multtping  Scnar  Decisions 

1.  INTRODUCTION 

This  report  is  a discussion  of  the  average  total  cost  of  making  a "target  present"  or 
"target  absent"  decision  from  the  information  contained  in  a sequence  of  sonar  pings.  This 
cost  is  useful  for  comparing  competing  decision  algorithms.  The  total  cost  is  computed 
as  a sum  c ti)  the  costs  of  terminal  decisions,  and  (2)  the  cost  of  pinging. 

The  d:  icussion  is  divided  into  two  parts . The  first  part  is  the  derivation  of  a formula 
for  the  total  cost  of  using  any  sequential  decision  algorithm  in  a multiping  sonar  receiver. 
The  derivation  presented  here  is  a straightforward  extension  of  existing  cost  computing 
procedures  for  one  possible  target  track  to  the  case  of  many  possible  target  tracks.  The 
second  part  is  a description  oi  three  sequential  decision  algorithms,  and  a comparison  of 
their  merits  based  on  costs,  error  probabilities,  and  ease  of  implementation. 

2.  COMPUTATION  OF  AVERAGE  TOTAL  COSTS 

The  multiping  sonar  receiver  makes  its  "target  present"  or  "target  absent"  decisions 
on  the  basis  of  the  returns  from  a sequence  of  sonar  pings.  We  assume  a receiver  which 
examines  each  possible  target  track  after  every  ping,  and  chooses  between  three  alternatives 
for  each  track:  (1)  decide  "target  present"  in  the  track,  (2)  decide  "target  absent"  in  tho 
track,  or  (3)  make  no  decision  and  ping  again. 


To  evaluate  this  receiver,  we  consider  the  following  experiment.  The  sonar  system 
continues  to  ping  until  a decision  of  "target  present"  or  "target  absent"  is  made  for  ever* 
track.  The  receiver  is  then  evaluated  on  the  basis  of  the  average  total  cost  of  performing 
this  experiment. 


in  order  to  derive  a formula  for  this  average  total  cost,  we  introduce  tM  following 
fjoUiS  ion . 

k = the  track  number 

k „ = the  total  number  of  tracks 

n t n v 


a random  variable  which  takes  or  the  value  1 when  the  decision  Is 
"target  present"  in  track  k;  and  0 when  the  decision  is  "targe?  absent*1 
In  track  k 
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a random  variable  which  taken  on  the  vaiue  1 when  e targat  is  actually 
present  in  track  k;  and  0 v'hen  a target  is  actually  ox&ent  in  track  k. 

the  cost  associated  with  the  experimental  outcome  i^ 

the  fixed  cost  of  a ping,  and  is  assumed  the  same  for  each  ping  (this 
cost  would  include  things  like  the  cost  of  alerting  the  target  and  the  energy 
required  to  generate  the  ping) 

the  cost  of  processing  a track  for  each  ping,  and  is  assumed  the  same 
for  each  track  and  each  ping  (unlike  c^,  this  cost  disappears  if  a decision 
is  made  which  terminates  a track) 

a random  variable  whose  value  is  the  total  cost  of  the  experiment 
the  a priori  probability  of  the  event  ij^ 

a random  variable  whose  value  is  the  number  of  pings  required  to 
make  a decision  in  track  k 

*i  random  variable  whose  value  Is  the  number  of  pings  required  to 
perform  the  experiment,  i.  e. , to  make  a decision  In  ail  tracks  (n  is  the 
largest  value  of  nk  at  the  end  of  the  experiment) 


The  average  total  cost  of  the  experiment  is  the  sum  over  all  tracks  of  the  average 
decision  cost  for  each  track,  plus  the  cost  of  ping’^g  required  to  perform  the  experiment. 
The  average  decision  cost  is 

kmax  1 1 

<decision  cost>  * E EE'WW  " 

k«*l  j=0  1=0 


The  average  pinging  cost  is 

k 

k „ 

max 

max 

<plnging  co8t>  = <Cjn  + 

c \>  “ cf  <n>  + S 

k*l 

k=l 

(2) 


The  average  total  cost  of  the  experiment  is  therefore  given  by 


max 

«>  ■ £ 
k=l 


1 1 


E £ «<w  »<w + cP  <v 

V*  ‘k'° 


max 

£ 

k»i 


Ml" 


<Cjt>  + Cj  <n>  (3) 


where 


Cjt>  * the  average  track  cost 
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This  result  shows  that  we  must  consider  the  average  track  cost  <Cjt> , plus  the  total 
cost  of  pinging,  when  comparing  the  cost  of  various  multiping  algorithms. 

The  expression  for  the  track  cost  can  be  further  broken  down  by  noting  that 
<nk> 

"k 

=EEnkpMk]v 

p«k> 

j 

Therefore 


<ck> 


E 

V0 


EcdWp[w+cPw 

Hr0 


p(U 


(5) 


In  order  to  study  the  average  track  coat  for  a single  track,  we  drop  the  subscript  k 
and  define  the  following. 


p [i^  = l/]k  = 0]  = probability  of  a false  alarm  = a 


p[ijt  = 0/j|t  = l]  = probability  cf  a false  dismissal  = 0 


pft  — l/j  -1]  = probability  of  a true  alarm  (detection)  = 1 - a 

k k 


p [tk=sO/jk=o]  “ probability  of  a true  dismissal  * 1-0 

p (||^  = 1)  = a priori  probability  of  a target  present  = P 

p(jk  = 0)  = a priori  probability  of  a target  absent  = q 


cd(l0)  > cd(00) 

cd<01)  > cd(ll)  (7) 

Therefore  all  the  terms  in  Equation  (6)  are  nonnegative . 

3.  MULTIPING  SONAR  SEQUENTIAL  DETECTION  ALGORITHMS 

In  thiB  section  we  describe  three  sequential  detection  algorithms,  and  compare  their 
relative  merits  on  the  basis  of  costs,  error  probabilities,  and  ease  of  implementation.  The 
throe  algorithms  are  optimum  sequential,  Wald  sequential,  and  truncated  Wald  sequential. 

Each  of  these  algorithms  are  developed  In  the  literature  considering  only  one  track. 
Equation  (3)  shows  how  to  evaluate  these  algorithms  in  the  multiple  track  situation. 


iwwsb 


•$r£3£l 


a.  Optim/m  Sequential 

The  optimum  sequential  test  is  described  by  Birdsall  and  Roberta  (Ref.  2).  This 
test  assumes  that  a maximum  number  of  pings  »mftX  are  allowed.  Within  this  constraint, 
the  algorithm  minimizes  the  average  track  cost  . In  the  multiping  sonar  case  there  are 
usually  a large  number  of  tracks  and  a small  maximum  number  of  pings.  This  combination 
makes  it  virtually  certain  that  pings  will  be  required  to  make  a decision  In  evory 
track.  Thus  <n>  = n , and  the  average  total  cost  given  by  Equation  (7)  will  be  mini- 

Dfl&X 

mired  if  the  average  track  cost  is  minimized,  i.  e. , if  the  optimum  sequential  test  Is 

. S'? 

The  optimum  sequential  test  compares  the  likelihood  ratio 


uBe<l  in  each  track. 


probability  of  the  received  sonar  echo  given  a target  present  ^ 
probability  of  the  received  sonar  echo  given  a target  absent 


& 


(8) 


to  two  thresholds  after  each  ping.  If  X is  above  the  upper  threshold,  the  decision  "target 
present"  Is  made,  and  the  test  terminates.  If  X is  below  the  lower  threshold,  the  decision 
"target  absent"  is  made,  and  the  test  terminates.  If  X is  between  the  thresholds,  no  deci- 
sion Is  made,  and  the  test  continues  to  the  next  ping.  The  thresholds  move  closer  to  each 
other  on  each  successive  ping,  until  they  coincide  on  ping  nmax*  which  forces  a declatcn  t'r.d 
terminates  the  test.  At  each  stage  of  the  test,  the  tradeoff  Ib  between  taking  another  ping  to 
Improve  error  performance  and  «.he  cost  of  pinging  again  to  obtain  such  an  Improvement. 

The  advantages  of  the  optimum  sequential  test  are:  (1)  the  ability  to  set  an  upper 
limit  on  the  number  ol  pings,  and  (2)  the  fact  that  costs  are  less  than  any  other  sequential 
test  with  this  feature.  The  disadvantages  of  the  optimum  sequential  teat  are:  (1)  the  diffi- 
culty of  computing  the  thresholds,  which  are  based  <xt  uncertain  costs  and  prior  probabilities, 
and  (2)  the  difficulty  in  relating  threshold  settings  to  posterior  error  probabilities  a and  8, 
which  are  often  the  only  meaningful  performance  measures  when  costs  and  prior  probabilities 
are  unknown. 

b.  Wald  Sequential 

The  Wald  sequential  teat  (Ref.  3)  is  the  special  case  of  the  optimum  sequential  test 

for  which  n is  infinite,  i.e.,  there  is  no  limit  to  the  number  of  pings  allowed.  Com- 
m 8X 

pared  to  any  other  test  with  the  same  posterior  error  probabilities  a and  JS,  the  Wald  Se- 
quential test  minimizes  the1  average  jwmbo*  of  pings  <»(0>  » and  <»(!>  > required  to  make  a 
decision  in  each  track  (Ref.  4).  TbU  minimisation  Is  lndeponsksl  of  the  assigned  prior 
probabilities  p and  q,  and  the  ©opts  e5  and  c . Thus  for  specified  a and  $,  the  Wald  sequen- 
tial test  minimizes  sverage  track  coots.  However,  the  tack  of  an  upper  bo«r<d  on  tbs  mvmV  ;•  r 
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of  pings  could  be  a considerable  disadvantage  when  the  algorithm  is  being  used  for  a large 
numb  r of  tracks.  The  value  of  <n  > could  be  very  large,  and  the  cost  term  Cj  <n  > te 
Equa  tor  (3)  migh.  be  intolerable. 

In  the  Wald  sequential  test  the  likelihood  ratio  X given  in  Equation  (1)  is  compared 
with  two  fixed  thresholds,  and  decisions  made  in  the  same  manner  as  in  the  optimum  se- 
quential test.  These  thresholds  can  be  set  to  achieve  any  desired  posterior  error  proba- 
bilities a and  0,  or  they  can  be  set  to  minimize  average  costs. 

The  advantages  of  the  Wald  sequential  test  are:  (11  no  other  test  has  fewer  average 
pings  to  decision  per  track  for  the  same  posterior  error  probabilities  a and  ft,  Imd^Z^Ke 
thresholds  are  simply  related  to  a and  /3,  and  remain  constant  throughout  the  test.  The 
main  disadvantage  of  the  Wald  sequential  test  is  the  lack  of  an  upper  bound  on  the  number  of 
pings,  which  can  be  quite  costly  when  there  are  a large  number  of  tracks,  and  a high  fixed 
cost  of  pinging  c^. 

c.  Truncated  Wald  Sequential 

The  truncated  Wald  sequential  test  is  the  same  as  the  Wald  sequential  test  except 

that  an  upper  limit  n is  placed  on  the  number  of  pings.  If  a decision  is  not  made  by  ping 
max 

number  nmjtx,  then  a single  threshold  is  set  and  a decision  made.  The  truncated  Waid  se- 
quential test  is  an  attempt  to  combine  the  simplicity  of  the  Wald  sequential  test  with  the 
limited  number  of  pings  feature  of  the  optimum  test.  If  the  maximum  allowed  number  of 
pings  is  large  enough,  it  is  reasonable  to  expect  that  the  performance  of  the  truncated  Wald 
sequential  test  will  be  approximately  the  same  as  the  Wald  sequential  teat. 

Comparing  the  truncated  Wald  sequential  test  with  the  optimum  sequential  test  and 
the  Wald  sequential  test,  the  truncated  Wald  sequential  test  (1)  has  a higher  average  track 
coat  than  either  of  the  other  tests,  (2)  has  a higher  average  total  cost  than  the  optimum  se- 
quential test  with  the  same  r.  , (3)  can  be  made  to  have  a lower  average  total  cost  than 

ill  8X 

the  Wait  sequential  test  In  most  situations  where  c^  is  much  higher  than  c , and  (4)  is  simp- 
ler to  implement  than  the  optimum  sequential  test.— 
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Technical  Memo  No.  2 
January  2,  1970 
D.  L.  Jordan 

Contract  No. : N00024-69-C-J.268 

TUI"!  Kalman  Filtering  Applied  to  Track  Detection 

In  any  track  detection  function,  the  problems  of  detection  and  estimation  must  be 
addressed.  The  best  estimate  of  a potential  track’s  future  position  and  velocity  are  needed 
for  association  of  die  many  target-Uke  returns  with  their  individual  tracks  as  they  are  ob- 
served. This  technique  requires  placing  a gate  (search  area)  around  the  predicted  position 
of  a particular  target,  and  if  a signal  return  falls  within  the  gate,  it  Is  associated  with  that 
taixet.  The  optimum  size  and  shape  of  the  gates  require  a solution  to  the  *>  called 
association  probier  i . we  have  used  the  gain  matrix  from  the  Kalman  filter  along  with  a 
velocity  sensitive  component  tu  aet  the  gate  size.  The  gate  size  versus  ping  number  with 
and  without  this  velocity  component  appears  in  Figure  1 . 

The  predicted  position  and  velocity  is  needed  for  both  the  positioning  of  the  gate  and 
the  detection  of  the  target.  In  this  technical  memo  the  optimal  processing  of  the  sonar 
returns  to  extract  range  and  velocity  information  will  be  reviewed  using  the  approach 
first  described  by  Kalman  (Refs.  1,2). 

Tho  Kalman  filter  will  be  applied  to  a linear  system.  A linear  system  ir  one  that  obeys 
the  principle  of  superposition.  Physically,  superposition  implies  tLat  the  presence  of  one 
excitation  does  not  affect  the  responses  due  to  other  excitations;  there  are  no  interactions 
among  responses  of  different  excitations  within  a linear  system.  In  the  analysis  of  the 
combined  returns  of  different  excitations,  one  can  begin  by  considering  each  individual 
cause  and  effect  separately  as  though  all  other  causes  were  absent  and  then  sum  over  all  the 
individual  excitations  and  corresponding  effects.  linearity  as  applied  to  the  random 
variables  y,  y (measured  range  and  rate  rate)  means  that  as  Independent  measurements 
they  will  define  the  total  measured  state  oi  the  system.  Wife  this  in,  mii.d  the  set  of  linear 
differential  equations  describing  the  p’  ysical  system  can  be  written  as 

dx. 

“j-  = f u>  P> 


Figure  1.  Gat/;  Sizes  Used  in  Correlating  Plng-to-Ptng  Returns 


| f where  x.  is  the  set  of  random  vectors  that  d ‘flue  the  states  of  the  system  and  u is  related 

f ’[  to  some  random  noise  generator  or  forcing  run-  . ion  that  adds  a noise  component  to  the 

I j measured  states.  The  measured  state  vector  is  then  expressed  as, 

5 r 

| j yt  = hxa  + v 


The  three-state  case  Is  represented  by 


In  matrix  representation  these  equations  become 


where  x.{  is  the  target’s  acceleration.  For  the  three -state  case,  the  state  vector  can  be 
shown  to  have  the  following  form. 


Having  described  the  equations  of  motion  for  the  linear  dynamic  system,  some  pro- 
perties of  the  random  variables  that  are  required  by  the  filter  equation  will  now  be  given. 

As  mentioned  above,  while  defining  the  random  variables  in  Equation  (2),  u{t)  Is  a 
gnusstirsn  random  vector.  This  means  that  the  sequence  of  random  vectors  x(t  - 1), 
xft  i 1).  generated  by  Equation  (2),  is  it  Gauss-Markov  sequence.  Therefore,  the  best 
estimate  of  a future  state  can  be  made  without  knowledge  of  all  previous  history.  Since 
the  development  of  the  equations  in  the  Kalman  filter  depended  or  this  property,  it  is  sot 
surprising  that  this  filter  gives  a best  estimate  of  a future  state  based  on  the  last  obser- 
vati.m  of  the  states.  To  complete  the  picture  of  the  model  for  the  system,  since  «(t)  is 
gauss inn.  the  sequence  of  random  vectors  uft  - 1),  v.jt),  and  u(t  + 1)  are  normally  distri- 
buted such  that  the  cross -covariance  matrix.  . 

eov  [ e(tj),  u(f,2)]  - 0 

This  Implies  tha*  u(t.)  »nd'#(l,,)  are  Independent. 


- 


- 


Furthermore,  It  is  assumed  that  £[u(t)]  = 0.  Since  this  noise  generator  is  present, 
the  output  vector  y(t)  of  the  measured  signal  alone  must  be  modified  by  some  additive  white 
noise (gaussian),  v.  with  E[v(t>]  = 0,  Efv^)  ’(tg)]  = 0,  and  E[v(t)  v(t)T]  = covfrftj)]  * R(t). 
Physically,  this  means  that  v(tj)  and  v(t2)  are  independent  and  that  a considerable 
amount  of  inaccuracy  is  present  in  the  measurement  of  the  signal.  In  fact,  the  further  one 
looks  back  into  the  state's  history,  the  larger  the  indeterminacy. 


The  filtering  problem  then  is  defined  as  giving  the  best  estimate  of  all  states  of  the 
syntem  based  on  the  noisy  observation  of  the  observable  states.  Since  the  system  is  linear, 
the  estimated  states  are  expressed  as: 


In  + K*n-Hin> 


(3) 


where  K(t)  is  a continuously  updated  weighting  factor  imposed  on  the  variance  between  the 
observed  and  predicted  states  and  as  yet  has  not  been  defined  explicitly.  *„(*)  the  best 
estimate  of  x^ft)  based  on  the  current  observation  yn<t). 

In  order  to  determine  K(t)  it  is  necessary  to  minimize  the  mean  squared  estimation 
error 

L - E[(x  -£>  • (x  — x ) r ] (4) 


Since  the  errors  are  gtuisslan,  minimizing  Equation  (4)  will  produce  an  optimal  weight- 
ing (gain)  matrix,  it  is  in  this  sense  that  the  Kalman  filter  is  optimal;  however,  this  pro- 
perty alone  will  not  necessarily  achieve  considerably  higher  levels  of  performance  than 
some  less  "optimal"  Altering  method.  By  performance  other  authors  usually  refer  to  the 
rate  at  which  the  elements  of  the  covariance  matrix  decrease.  In  other  words,  with  in- 
creasing time  the  estimates  of  the  system  are  said  to  improve.  For  the  case  when  there 
are  several  data  sample  sets  however,  the  updated  estimate  of  the  gate  sizes  ("updated" 
here  means  that  the  size  of  the  gates  is  decreasing)  can  in  effect  be  reducing  fee  probability 
of  declaring  a false  track,  by  eliminating  certain  nonconsistent  tracks.  At  the  same  time, 
if  the  target  motion  is  not  consistent  w'fh  the  system  equations  of  motion,  the  detection 
probability  of  the  system  will  be  decreased. 

Using  Equation  (3),  the  loss  function  becomes 

E{(x  -x)  • (x  - x)1)  - E[(x  -x)-  (x  -x)T]  + E[(K<y  - Hx))(K(y  -Hx))TJ 

“ECKfr  -Hx)(x  -x)TJ  - Efcx  - x)  KT(y  - H x)T] 
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Using  the  following  equations,  y = Hx  + v,  and  noting  that  E[(x  - x)  v ] = 0,  the  lose 
function  now  becomes 

E(<x  - x)  (x  - x)T]  = P - KHP  - PHT  KT  + K(HPHT  ♦ R)  KT 

In  order  to  minimise  the  loss  function  we  set  the  derivative  of  L,  w • r * t • K,  dL/UK, 
equal  to  zero.  In  so  doing, 

-2PHT  * 2(HPHT  + R)  K = 0 

Solving  this  equation  for  the  gain  matrix, 

T T -1 

K = PH  (HPH  + R) 

With  this  value  for  the  gain  matrix  the  estimated  covariance  matrix  becomes 

A 

P - <!  - KH)P 

In  order  to  complete  the  Altering  problem  a recursive  relationship  for  the  conditional 
covariance  matrix  P(t  + 1,  t)  must  be  derived. 


P(t*l,  t)  «E[(x<t  + l) + 0),  (x<t  + l) -x<t+ 1,  t»)TJ  (6) 

■ but 

| *<H  1>  -•!»  ♦ t.  t)  x(t)  + Aft  + I,  t)u(t)  <7) 

1 

and 

+ *»-*«» + i. ict  ♦ 1. 1)  (8) 

■ 

| After  substituting  Equations  (7)  and  (8)  into  Equation  (6),  the  predicted  covariance 

matrix  becomes 


P(t  + l)  - *(t  + 1,  t)  P(t)  *(t  + 1,  t)T  + AE  [u  • uT]  AT 

For  the  two-state  case  where  the  model  error  is  assumed  constant,  the  predicted  co- 
variance matrix  fee-cornea 


With  these  values,  the  following  can  be  defined. 


1 


n 


*v<t) 


22 


<V-r  " P12 


21 


The  gate  slaes  (see  Appendix  B>  for  this  model  appear  in  Figure  1.  Curve  A contains 

2 2 
bo  model  error,  !.*.,  <ry  =0.  Curve  B contains  the  model  error  contribution,  «r  , but  its 

value  is  obtained  h am  the  predicted  range  rate 

For  the  two-state  case  with  nonconstant  model  errors,  the  predicted  covariance  matrix 
2 

has  the  form  (set  it  = 0). 


P(t  ♦ 1)  - 


» P 
r 11  ^12 


P21  P22 


♦ 


where  x^  is  the  expected  range  rate. 

Equation  (2),  as  it  ataitda,  does  not  give  tie  predicted  estimate  of  the  system.  Since  on 
the  average,  expected  values  will  result  the  predicted  estimate  of  the  state  vector  is  given 

by 


X(t  4 1)  » *(t  + 1,  t)  X (t)  + A(t  ♦ ;,  t)  E [u(t>] 


However,  since  E[u(t)]  > 0 for  all  t,  than 

~ A 

X(X  + l)  *>  *(t  + 1,  t)  jr(t) 

A summary  of  t»«s  filter  equations  is  given  in  Figure  2. 

In  order  to  start  fJse  filtering  problem  some  Initialisation  of  the  covarlroco  matrix, 
P(t<t),  for  t-!c  best  estimates  of  the  initial  states  of  the  model  must  be  made,  t.e.,  some 
systematic  selection  of  the  Initial  values  of  the  covariance  matrix.  The  satecfiiosi  of  these 
elements  depends  open  a knowledge  of  the  problem,  t.e.,  the  track  dynamics,  Unfortun- 
ately, a Hostile  target  doss  not  cooperate  l»  revealing  Us  maneuvers.  In  this  Instance,  It 
seems  that  the  heat  selection  would  fee  to  set  the  off  diagonal  elements  to  aero.  This  is 
reasonable  because  knowledge  of  the  Initial  range  («n  either  fee  first  ping'  or  first  return) 
wifi  in  no  way  provide  tefomaSoa  about  the  initial  velocity,  i.e.,  Initial  rang*  And  voice  fly 
are  tmeombted.  The  mate  diagonal  terms  arc  act  so  that  the  predicted  covariance  will 
set  a gate  about  a fMroot  target. 

v n 

3k  / 


It  should  he  noted  that  the  covariance  matrices  are  predetermined,  t.e.,  their  values 
are  independent  of  any  measurements.  Clearly  this  means  that  the  gates  are  updated 
^narrowed)  for  both  false  and  fading  targets.  In  order  to  rectify  this  insensitivity  to  the 
roal  world,  as  a preliminary  improvement  to  the  covariance  matrix,  the  calculated  pre- 
dicted range-rate  is  used  la  the  matrix.  Qualitatively,  this  seems  to  have  reduced  the 
root -mean  squared  error  between  the  predicted  and  measured  state  vector  more  quickly 
than  when  no  velocity  terms  were  used  in  the  model  errors.  More  work  on  this  point  is 
presently  in  progress. 

In  conclusion,  one  of  the  main  functions  of  the  Kalman  filter  is  to  update  the  measured 
state  vectors  obtained  from  previous  instants  of  time  to  a set  of  values  corresponding  to 
the  present  instant.  These  updated  measurements  are  used  along  with  the  current  meas- 
urements to  form  an  optimum  estimate  of  the  state  vector.  Since  a random  disturbance 
term  exists  in  the  canonical  equations  of  motion  there  Is  an  uncertainty  in  updating  the 
measurements  of  the  previous  instants.  The  Kalman  filter  takes  this  fact  into  account  by 
assigning  weighting  factors  to  the  contributions  of  the  previous  Instants.  The  fact  that  ( a 
mean  squared  estimation  error  is  minimised  and  that  the  Kalman  filter  uses  matrix  weVf'- 
Ing  factors  so  feat  each  component  in  the  state  vector  can  be  weighted  individually,  acco'.jv.a 
in  part  for  its  superiority  as  a filter. 


APPENDIX  A 


DERIVATION  OF  THE  TRANSITION  AND  DISTRIBUTION  MATRICES 
FOR  THE  TWO-  AND  THREE-STATE  CASES 


The  two-state  ease: 


x = x„ 


x2  ” Ux 


Integrating  Equation  (10)  first  w ■ r • t • "t" 
x2(t  +1)  = ux  T ♦ Cj 

and  Equation  (11)  becomes 
x2(t  + 1)  = ux  T + x2(t) 


Substituting  Equation  (12)  into  Equation  (9)  and  integrating  over  t,  keeping  Xg(t)  constant. 


Xj(t  + 1 


“x  ^ 


+ x2(i)T  + C2 


Again,  applying  the  boundary  conditions, 
u T 

*j(t  ♦ 1)  - + ^(t)  T + Xj(t)  (H 

where  Xj(t)  is  some  initial  position  at  T * 0.  Therefore,  the  canonical  equations  of  motion 
become 


Xjtt  + 1)  = x + T^  + 


x_(t  + !)  = 0 4 x.  + u.  T 

In  order  to  put  these  equations  into  matrix  natation,  first  examine  the  case  where  = 0, 
i.e.,  «he re  are  no  meded  errors. 


Xj(t  4 1)  - Xj  + Tx 


x2(t  4 1)  -*  0 4 x2 
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Clearly,  these  equations  in  matrix  notation  become 


il7) 


From  the  comparison  of  Equation  (15)  with  Equation  (16),  It  can  be  seen  that  another 
matrix  must  be  added  to  account  for  the  model  errors.  These  can  be  represented  by  the 
following  matrices. 


Model  Contribution  = 


0 


= A(t  + 1,  t)  u(t) 


(IS) 


The  derivation  of  the  equations  in  the  three -state  case  follows  the  exact  same  procedure. 
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